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Assuming that gluon condensates are important even above the deconfining phase transition, we 
develop a model for the dilepton yield from a quark gluon plasma. Using a simple fire ball description 
of a heavy ion collision, and various estimates of the strengths of the gluon condensates, we compare 
our predicted dilepton yields with those observed in the CERES and HELIOS experiments at CERN. 
The simple model gives an adequate description of the data, and in particular it explains the observed 
considerable enhancement of the yield in the low mass region. 

I. INTRODUCTION AND SUMMARY 

Extensive lattice simulations have shown that QCD exhibits a phase transition at a temperature of the order of 
the infrared scale A. This transition is believed to be a confinement-deconfinement transition for heavy quarks, 
and a chiral transition for light quarks. The high temperature phase is characterized by a small coupling constant 
g(T). Unfortunately, perturbation theory in g(T) is plagued with infrared singularities since the magnetic sector is 
non-perturbative at all temperatures. 

High temperature lattice jl]|2| and instanton simulations || seem to concur on the fact that in the vicinity of 
the critical temperature, T c , the chromoelectric and chromomagnetic condensates are still substantial. A number of 
theoretical arguments, both perturbative [Q and nonperturbative ||, suggests that in the Euclidean formulation, the 
high temperature phase of QCD exhibits strong A4A4 correlations. It has even been argued that the global color 
symmetry might be spontaneously broken with the appearance of an A4 condensate ||||. Also, lattice simulations 
show that chiral symmetry is restored at the critical temperature, so that the quark condensate can be set to zero 
above T c . 

In the present analysis we will incorporate these three important observations that are pertinent to QCD near and 
above T c ; the quark condensate vanishes, the chromoelectric and chromomagnetic correlations remain non-zero, and 
strong correlations develop in the static part of the A4 gauge field that can be described by a local A\ condensate. In 
the spirit of the QCD sum rule method, we will not address the complicated dynamical questions of how the gluon 
condensates come about, but simply assume that they exist. We will also restrict our discussion to only those with 
dimension 2, {^-Af), and 4, (^-E 2 ) and (^-B 2 ), following an earlier suggestion by two of us |J. The value of the 
condensates can only be estimated using non-perturbative calculations such as lattice QCD, and at high temperature, 
the scale is expected to be set by the mass gap in the magnetic sector which is ~ g 2 T. 

That chromoelectric and chromomagnetic condensates could be present in the high temperature phase is not very 
surprising, and their size can be roughly estimated using the bulk energy and momentum measurements at finite 
temperature on the lattice. The analysis in g] gives the values : {^-E 2 ) ~ (^B 2 ) ~ (200 MeV) 4 

The presence of an A 2 condensate is more controversial Q , so we will give a summary of the arguments presented 
in M for why such a condensate should be considered. The basic observation is that the Wilson line (or Polyakov loop) 
operator W — TrPexp(i dx^ ^4(2:4)), is a gauge invariant operator that is a good order parameter for the pure glue 
theory. Although this operator is non-local, it is possible to relate it to a string of local operators by choosing the 
static gauge d^A^ = 0. For example, in an SU(2) theory we have (W) = (1 + cos^A^)), and one could now try to use 
factorization to extract (A 2 ). Here we will simply equate the condensate to the magnetic mass, i.e. we will typically 
use, {^A\)/T 2 ~ ^-(g 2 ) 2 ~ 16na^, which is about 0.4 for a s ~ 0.2 This estimate is also supported by a lattice 
calculation in pure Yang-Mills Indeed, an analysis of the gluon- propagator in four-dimensions in Landau-gauge 
gives a 2 ((gA/C) 2 ) ~ 1/3 on a 4 x 8 3 lattice, which is equivalent to (^A 2 )/T 2 ~ 0.4. A direct reconstruction of this 
condensate from the dimensionally reduced theory has proven to be more subtle ]§[ . 

In this letter, we present a calculation that suggests that the gluon condensates strongly affect the emission rates 
of low mass dileptons from the plasma. Our calculation will be very naive. We will simply assume that the quark- 
antiquark annihilation processes that give rise to the dileptons occur in the presence of low momentum fluctuating 
background gluon fields. By averaging over the background fields, we get leading corrections to the emission rates 
proportional to the various condensates. Note that in doing so we take the concept of a fluctuating background field 
at face value and perform the calculations directly in the time-like region assuming that the temperature is high 
enough for perturbation theory to be useful. Thus our method is very different from the QCD sum rules, where the 
background field method is just a convenient tool to calculate coefficient functions in the operator product expansion 
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which is then used in the deep Euclidean region.^] However, since the emission rates can be expressed in terms of 
current-current correlation functions we can use the standard QCD sum rule techniques to calculate the condensate 
contributions. How this is done in detail, especially when it comes to gauge choice, is explained in fij. To calculate 
the contribution from the A 2 condensate we will use the following observation: A constant potential A4 is the same 
as an imaginary (color) chemical potential, fi a , and thus one can obtain the term ~ (A 2 ) by using a grand canonical 
ensemble for the quarks with a chemical potential /i a = iA%, averaged over the arbitrary directions, a, in color space, 
and pick up the term ~ /i 2 . (We have also checked that this procedure gives the same result as the standard fixed 
point gauge methods used in 0.) 

At this point we want to stress that although our approach is inspired by the above theoretical arguments, it 
is basically phenomenological, and should primarily be judged by how successfully it can accommodate data. The 
theoretical uncertainties are not limited to the assumptions we made about the non-perturbative corrections to 
correlations functions in a high temperature quark gluon plasma. We also have to make several assumptions about 
the physics of the heavy ion collision, and we will use an oversimplified model, namely that of an expanding fireball 
containing only the plasma phase. The main motivation for using such a simple model (the details are given in section 
III below) is to demonstrate that the gluon condensate effects are important. Hopefully more realistic models e.g. 
including regions with mixed or hadronic phases, would also exhibit this sensitivity to the condensates. Also in the 
transition regime, the pure and mixed quark condensates may still be relevant. 

Having said all this, we will demonstrate that a naive fire ball model combined with a simple parametrization of the 
production amplitudes based on the above ideas, and with the condensate parameters at the values suggested by the 
scale arguments, is quite successful in explaining the dilepton data from the CERES and HELIOS relativistic heavy 
ion experiments. The rest of this paper, where this claim is substantiated, is organized as follows: In section 0, we 
calculate the effe cts of the leading order condensates on the dilepton emission rates from a thermalized quark-gluon 
phase. In section III, we use a simple fire-ball model to estimate the dilepton yield, and compare to the CERES and 
HELIOS experiments. In section IV, we show that our emission rates are also consistent with the recent analysis of 
the CERES results using different p t cuts. 



II. EMISSION RATES 



In a translationally invariant, but non-perturbative, quark-gluon phase, the rate R of dileptons produced in a unit 
four volume is directly related to the electromagnetic current-current correlation function ]Io| , [rT[ | . For a pair of leptons 
with mass mi and momenta ^1,2) the rate per unit invariant momentum q = pi + P2 is , 

dR a 2 / 2m? \ / 4m 2 x 1/2 



d*q ^W\ l + Z f)V-^F) W(q) ' (1) 

where a = e 2 /An is the fine structure constant, and 

W(q) = J d 4 xe~ tq x TT (Pe iA S Q V (H - F)/7 >(a:)J M (0)) . (2) 

Here eJ M is the hadronic part of the electromagnetic current, H the hadronic Hamiltonian, F the free energy, T the 
temperature, and Q a a global color generator. (As mentioned above the direction, a, is arbitrary but fixed.) W(q) is 
related to the retarded current-current correlation function by, 

W(q) = eq0/ l _ i lmW R (q) (3) 



where 



W fl (g) =i /d 4 a;e i9X Tr (Pe iA 4 Q %-( H - F )/ T R*(J^(a;)J^(0))) , (4) 



1 A work similar in spirit to ours is the calculation by Patkos and Sakai |12| of the nonperturbative magnetic contribution to 
the e + e~ annihilation rate. That calculation was inspired by the so called Copenhagen vacuum based on fluctuating magnetic 
fluxtubes. 
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and R* is the retarded product. The leading contribution to the retarded function is condensate independent, and is 
just the "Born" qq annihilation term, which for massless quarks becomes, 



MW * = ir c £ 



q 

q— u,d,s 




2T fn+ 



(5) 



where N c is the number of colors, e q the charge of the quarks, and n± = \/ (e^ qo± ^/ 2T + 1). At T = this is just 
the i?>3 ratio well known from e + e~ annihilation 

^>3(g 2 ) = E i , = ImW oVo = ■ (6) 

i 9 

Using only the three light flavors, i?>3 is 2. 

As mentioned above, the contribution of the A4A4 condensate to the retarded function can easily be calculated 
by picking up the term quadratic in the imaginary colored chemical potential A\. The result is 



ImW* = ± Uj2ei\ q 2 [(^A 2 ) (n + (l - n + ) -„_(!- 



2 ■ 

(7) 



According to the estimates in the previous section, (5 = (^Al)/T 2 is of the order of 0.4. Note that (^) has the same 
sign as (|5|) and will thus increase the dilepton yield. In the next section, we will see that this effect can be rather 
dramatic in the invariant mass region studied in the CERES experiments. 

The contribution of the chromoelectric and chromomagnetic condensates to the retarded function can be calculated 
using the methods in |Q, and the result for q 2 > is, 



-wg -i UxA [(_i<**) + i<**» (±9 ( 



(8) 



The temperature dependence of (||) is the same as that of ([!]) and can be understood from the condition of detailed 
balance. Note that the emission rate is enhanced by a chromomagnetic field and depleted by a chromoelectric field 
as expected^], and also that all condensate contributions vanish at T = 0. 

III. RESULTS 

In Fig. [I] we show the dilepton emission rates at temperature T = 150 MeV, using the parameters, S = 
(^E 2 }/(200 MeV) 4 = (^B 2 }/(200 MeV) 4 j|. The contribution from the non-perturbative plasma, i.e. perturba- 
tive plus condensate contributions (dotted and dash-dotted lines) should be compared to the purely perturbative part 
(thick solid line) and, as a reference, to that of a gas of mesons (thin solid line) jl3| . The condensate contribution 
to the rates are sizable in the 500 MeV region, and the effect of the (A 2 ) condensate, /?, dominates that of the 
chromoelectric and chromomagnetic ones, S, for masses above 250 MeV. The dependence on both condensates /3 and 
5 is negligible at energies higher than 2 GeV, and above 2.5 GeV the perturbative plasma gives the same result as 
that of a hadron gas, as expected from duality. In the Dalitz region (< 200MeV) the effects of both the electric and 
magnetic condensates are appreciable and dominates the perturbative part by orders of magnitude. This observation 
maybe important for photon emission as well. 



2 Roughly speaking, the e + e pairs are repelled from each other in the CM frame in the presence of a chromoelectric field, 
and deflected in the same direction in the presence of a chromomagnetic field. 
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FIG. 1. Dielectron rates at T=150 MeV. = (^A 2 4 )/T 2 and S = £ 2 }/(200 MeV) 4 = {^B 2 )/(200 MeV) 4 . See text. 

In a heavy-ion collision the electromagnetic emission occurs from various stages of the collision process, and hence 
at various temperatures. In this part, we will focus on the emission rates in the low and intermediate mass region 
(up to 1.5 GeV). We will assume that the dileptons emanate from a simple fire ball composed of the nonperturbative 
plasma phase only. The fire ball expansion is assumed to be homogeneous, and we will use a standard parametrization 
of volume and temperature, 



V(t) = V 1 + 



t 

to, 

T(t) = (T l -T 00 ) e - t / T +T 00 . (9) 

The parameters in (|^) regarding the initial and freeze-out temperature will be taken from present transport calculations 
Jl4],[l5]]. We will not try to justify this point except by saying that it allows us to define a space-time volume for the 
"hot" spot in the heavy- ion collision that is consistent with the one used for an interacting meson-phase [p~3| — p~5|| - 
Hence to — 10 fm/c, T t — 170 MeV, = 110 MeV, r = 8 fm/c, and the value of Vq is absorbed into the over-all 
normalization constant -ZVqVo = 6.76 x 10~ 7 fm 3 . The freeze-out time will be set to £/. . = 10 fm/c. 

Using (^|) and the above rate, the final expression for the integrated emission rate per unit rapidity r\ and invariant 
mass M is 

dN/dndM , r f tf °- , Tr , s f d 3 q „, 2n dR 

' (M - N M / dtV{t) / ^A(q ,q 2 )— . (10) 

The acceptance function A enforces the detector cut pt > 200AfeV,2.1 < r\ < 2.65, and 9 ee > 35 mrad for CERES. 
In Fig. |^, we show the evolved emission rates for CERES. The enhancement is sizable and goes in the direction of the 
data in the 500 MeV range in CERES. Again, the nonperturbative QCD results are indicated by the dotted line for 
the condensates we quoted, the dashed-dotted line (without the E- and B-condensates) , while the perturbative QCD 
result is indicated by the solid line and the meson result by the thin solid line. 
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FIG. 2. Dielectron rates for CERES S-Au experiment. See text. 

For HELIOS-3, © can be used by integrating over n with the cut nit > 4(7 — 2^), nit > \J (2m M ) 2 + (15/ cosh(?7)) 2 . 
In Fig. ||, we show the evolved emission rates for HELIOS. For this energy range the effect of chromoelectric and 
chromomagnetic condensates are small irrespective of the specific value for the thermal condensates as indicated in 
Figs. |l| and 0. The enhancement is compatible with the HELIOS data, and slightly overshoots both the data and the 
meson gas in the 1 GeV region. 




1.0 1.5 2.0 2.5 

M (GeV) 



FIG. 3. Dimuon rates for HELIOS-3 S-W experiment. See text. 
IV. P T CUTS 

Another important check on the nonperturbative emission rates can be obtained by using the recent p t cut analysis 
by the CERES collaboration to probe the transverse momentum content of the dilepton emission rates in the fire- 
ball. Our results versus the data are shown in Figs. ^[^, with and without the chromoelectric and chromomagnetic 
condensates. The pt dependent effects caused by the thermal condensates are most dramatic in the low p t region 
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(Fig. ||), where there is a sizable condensate effect, mainly due to (3. The condensate contributions are in the right 
direction around the low mass dielectron enhancement, and give important additions to the perturbative result. 

In the high p t region, Fig. 0, the resulting condensate effects are rather small, and the curve is consistent with 
the present data. Therefore, the total transverse momentum content (Fig. ^) is quite accurately reproduced by the 
addition of the gluon condensates. 

In summary, we find it encouraging that not only the dilepton emission rates, but also the observed low mass-region 
pt enhancement is, at least partially, explained by our simple model, and we think that this further strengthens our 
proposal that nonperturbative condensate effects play a role in the plasma phase of heavy ion collisions. 
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FIG. 4. Dielectron rates for low p t of CERES Pb-Au experiment. 
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FIG. 5. Dielectron rates for high p t of CERES Pb-Au experiment. See text. 
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FIG. 6. Dielectron rates for the total p t of CERES Pb-Au experiment. 
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